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Destabilizing fluid forces on a whirling centrifugal impeller rotating in a 
volute have been observed (Ref.l). A quasisteady analysis neglecting shed vorticity 
(Ref. 2) or an unsteady analysis without a volute (Ref .3) does not predict the 
existence of such destabilizing fluid forces on a whirling impeller. 

The present report is intended to take into account the effects of a volute and 
the shed vorticity. We treat cases when an impeller with an infinite number of 
vanes rotates with a constant velocity Q and its center whirls with a constant 
eccentric radius e and a constant whirling velocity <u. 

Major assumptions are as follows: 

(1) The number of the vanes is so large that the impeller can be treated as 
an actuator impeller in which the flow is perfectly guided. 

(2) Flow is inviscid, incompressible and two-dimensional. 

(3) The eccentricity e is so small that unsteady components can be linear- 
ized. 

(4) Vorticity is transported on a prescribed mean flow, i.e., the operating 
point is near design flow rate. 

(5) The volute can be represented by a curved plate (see Fig.l). 
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SYMBOLS 


2 = Z+*y : 


z' = X'+tf' 
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nr : 

(UTr, ur 9 ) 

t : 

W : 


stationary frame with its origin 0 fixed to the center of the whir- 
ling motion 

translating frame with its origin O' fixed to the center of the 
impeller and its axes parallel to those of Z. 

inner and outer radius of the impeller 

vane angle measured from circumferential direction 

eccentricity 

angular velocity of whirling motion 

angular velocity of the impeller 

prerotation, Q: flow rate 

absolute velocity 

velocity relative to x'y* frame 

velocity relative to the impeller 

time, t=0 when 00' is in the direction of x 

circumferential mode number 


Subscripts 

1,2 • quantities at the inner and outer radius respectively 

Superscripts 

: steady component ~ ; unsteady component 


BASIC EQUATIONS 


Euler's equation in the rotating and translating frame fixed to the rotor is. 


dt* 




(!) 


where d/dt,_is the time derivative in the rotating frame and JJ^iwee w and 
0 X x^irQe 1 are the translational velocity due to whirling and the rotational 
velocity of the impeller. In an actuator impeller with an infinite number of vanes. 
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the effects of the vortlcity distribution on the vanes and the forces exerted by the 
vanes are represented by the vorticity V X v and the external force £ 
respectively in the above equation. For inviscid flow through actuator impellers, 
£ and w X (VX j) are normal to the vane surface and the component of Equation (1) 
parallel to the vane surface is 


dVj 

dt* 


+ 


a. / ur z Vt z 
M<S " 2 



- 0 


c> 


where 

Vs - hT ~ rJ2a>S(2 - ojS coi (&tfi-ujt) 

= WrCK , 0-9cn) r z / c r<s)*us»p(n) 

Vt z ~ f £*cj 2 i- zr(o£{2ccA(o-a>t) 

+ 12)/ *6 

Integrating Equation (2) along the vane surface we get the following total pressure 
increase. 

£t f &L * I m Ve'J*- H (£ r£*) Ifrz - CQ-0>)j*0*,(&+p -cot)d$ 

t Sto[( m t V# ) a>s(&-(Ot) r W/SM(8-i6V]f ( 3 ) 

where 

M * f * - - dr 

J, r Ac* per) 

Downstream of the impeller, where f=0, Equation (1) can be expressed as 

i / , V Vt 2 

%rx ( V* 2f) - 7 ^ f iVH t H -j 'j" "T" <*> 

By multiplying g*alz.-t gydy - l/l(efteS8 -&zsi*>9)d&- on both sides of the above equa- 
tion we get the following vorticity distribution at the outer radius of the 
impeller; 

^ I aHt I / d a \ / 

* “ MSrz * 6 

= f t Vr'c*b(&,-(t)£ ) - (&, -U)t) 

t ^cosCBroX) (Pi -cot) - ir,Cl<U#>(Bi -tot) j- 

f (x e*' M + Sco(&-tojfas(9rp -cot )dS 

' r Q-i&)Vo2 aj 


163 



ELEMENTARY FLOW COMPONENTS 


The flow tangency condition at the impeller outlet is. 


Viz = iz C2 ~ V /2 ait ^ (6 ) 

The first term on the right hand side of the above equation is cancelled by the flow 
component 

U~r ~ *• Vo? ~ ~2 Jl7~ { ® — tkfi ~~ (.7) 

and other disturbance components should satisfy the following equation 

ir&z * - v r i &>t/3z , m 

For the region outside the impeller we consider two types of elementary velocity 
disturbances satisfying Equation (8): 

First, consider the following velocity field: 


U ''~ A ' V ' ~ 27L I z'~ & + e ^ & £'~nv% (<f) 

This velocity field has a vortex T at z'=z'q and satisfies Equation (8) with no 

circulation/net flow around/from the impeller. Consider a vortex distribution 

Hs) = rs(s) + el^slsin wt + eI^(s)cos ut on the volute surface z (s)=z* +ee iwt . 

a a o o 

Assuming e << we get the following steady and unsteady velocity components 
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where z* is fixed to x’y'-plane. At z=z =z'+ee* wt * which is fixed to the volute, 

s 

we get the following expressions: 


u'-iirlis* 



i 

z s - 2,(S) 


rf e zifL 7 
z s ($?>«) - r/) J dS 


(12) 
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- ,jv . _ il TV r/e *&* , e i»t 2 ,(s>e* wt t B s e* at s .. 

U'-±V h s - 21L } 0 JsCS) Zt (Bs2o c t> - £*) 'ft + ft + r.* yi ^ 

In the same way the velocity component (7) has the following steady and unsteady 

i& s 

components at the volute surface z =r e “z'+ee 1 * 0 . 

s s 


Vr - l U»L = IJcT s { a ~ 2x1 K(nf 2 2 ® catfo )} 

*5 

W-i ir&h IW s i Q " 2rUrz(&n-5ffiC*tp2)}~e 

“s * 

Next consider the velocity field due to shed vorticity. 
vorticity is transported on a streamline 


If we assume that 


04) 

OO 

the 


w-ivS* (&-xT)/(nn 


we have the following elementary vorticity field 

5* = l<mCKizu>(t-$cr*-r?>) + *» (»-§*}Cr/ tk»} 

t £ s „ SM>l±toCt--$-(r*-a‘))i- /»(e--§-&}(r/ri)) } m) 

and the corresponding velocity field 

ITfW ~ A V~6W ~ (Teemin'/)* £$ 7 ? ?XV>)c0S (W81u>t> ) + 2 BO) 'Kco^xoo ) $A/A tW&itOt> ) O'?) 

where 


Zj ?m - Ren - i <6>fivi 


2lo> - ’Rim ~ l <®xo? 


Rp'n + ^ Riv * -j. [ F(r) + Grcr) ] , 0 e-n v @10 J = -F(r) t drir) 
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dcr) = j { j. (teo^crS- n 3 ) - m ( r*/n»j C r/ n )*"Ur 0 


Equation (17) includes the effects of the vorticity in r,<r<R. Since G(r) is 
convergent only for n£3 as R->®, we should artificially prescribe some 
appropriate finite value for R. Equation (17) does not satisfy Equation (8) and we 
should add a potential component 
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such that 


Ubn ~ = C Vr& ~ U&oi ) + C n ” * 1/gS ) 


av 


satisfies the boundary condition of Equation (8). We may now note that there are 
steady and unsteady components as follows: 

(i) Steady component; in Equations (16) -(19) we put u>=0 and represent the 
related quantities with superscript and suffix n (n=l,2...) for each 

mode * Hem etc * 

(ii) Unsteady component; correspondingly to the sign on u and the mode n, we 
represent the related quantities such that etc. 

On the volute z=z s =z'+ee* ft>t , the steady and unsteady velocity fields are expressed 
as follows: 


U'-iU'li , = ^ If + tomlttis] e~ 2&s 

■* m-i 

— 6 *• ( fan ~ ’**■ J feri./ (20) 

09 

U'-iV'k s /Z = 2- "a- C&StyQs) 

m-i 

+ (HsmZpm - Q*2i:*)('r>C£S'r>8 s ^ SimCfig -cot) e *** 


~ ■*" CimBrm ')COS l Y)&s + C ~Hcm (&s -(Ot,) Q 

~ (fit I r J 
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V> s l 
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In the region upstream of the impeller (r^r.,) , we consider a source Q and 
prerotation of strength F^ at the center of the volute. Then the velocity can be 
expressed in the series 


Vr -l rn'= [fiT WTse^ *vJ A/r> ^ S tAmcost ot)j 'j e LQr 


(22) 


— ■" ■" s s see c 

where A^=Ap r[ +iA TTi , +iA Tf>> A^A ^+ iA^ are complex constants. 

Now we have given all of the elementary flow components necessary for the 
construction of the entire flow field. Each of them contains several unknowns that 
are determined in the following sections. 
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BOUNDARY CONDITIONS ON THE VOLUTE 


The flow tangency condition on the volute surface is, 

SU) Cd>S (C<-COt) + Ir'CtXK* - U'/tlnOl -a V* = 0 (2$) 

where a is the angle between volute and x-azis. If we put 

h! - U' + UC CoS cot i- Us si * 1 col 

l/' - lf'-t ~Ul CoS cot + V$ sin cot 

we get the following conditions, 

ir'ama - R'&Imoi = o > (2<t) 

dZotSol - Ui/dlnoC - - SU> CeQCX , C2& 

ITs'cosot - Us - -£u)Oon<X. ( 26 ) 

The steady velocity u' ,v' is given as a sum of the velocity components in Equations 

(12), (14) and (20). The unsteady components u' , "v*' and ^u* , *u’ are the cosine 

c c s s 

and sine components of the velocity of Equations (13), (15) and (21). Equations 

(24)-(26) constitute integral equations for r s (s),I^(s) and I*^(s). 

CONTINUITY EQUATION 

The continuity equation across the impeller is 

tTrc n, 0 ,* =■— Vr cr z ,&) C21) 

where is the angle between corresponding leading and trailing edges of a vane. 
At the outer radius the total of the steady velocity components given by Equations 
( 7), (10) and (19) can be expressed as a Fourier series; namely, 

q . &? _ , 

ir/(n,&> -TiFr * 2 . C Vn» x>9 t cos nod ) 

* w-i 

_ w _ 

VZCr^e) - V$vsiM'r)& + lx e %ct>s'r7&) 
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and (19) can be 


In the same way, the unsteady components of Equations (11), 
expressed as 

lf/cr 2> &) = 2.2. i ( Vm cos m& t VmtoM'nd) cosuit + ( w9 r V^s^9>£>)s2 u>c]I 
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At the inner radius * =t * we expand Equation (22) in 0=0., rather than in 0 1 ; 

then. 
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V&cr h 9 ,=G+?o) = 22 S( VivcxwP + l%$><u»w 9 )c*svt+( 

From Equation (27) we get the following relations 


J 


vii, * (. ft/tf > V& 

C32) 

Vm * CtL/r,) v£> 


(33) 

Vm = CWr,) Vrm 

OO 

Vrm =cn/K) ir£ 


(3£> 

Vm= Ch/r,) v% 

06) 

&£ = (. fc/ r, ) Uri 


07) 

Equations (32) and (33) give the relations to determine A„ and 

tions (34)— (37) determine A^ q , A^ & , A® n and Aj^. n 

A In' 

and Equa- 

STRENGTH OF SHED VORTICITY 



If we use the expressions (28-31) 
components can be expressed as 

in Equation (5), the 

steady 

vorticity 

jr _ Tr c 

hem ~ q Ju'H ( £ Urn 

+ 
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where we have used because of the assumption on the transport of the 

vorticity. In the same way, if we express the vorticity by 


00 

'%(&>&) = 2 ^~ s 2 [c Ql c Casr)& r^smv>e) si*)U)ti-( C^cuSViQ + 7iB)a>Scot}} ( 40 ) 
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and use the expressions (28-31) in Equation (5), we can express £ in terms of the 
Fourier coefficients in Equations (28)-(31). Comparing Equation (40) and Equation 


we get the following 

relations. 



& -Sfr&t e?) 

(41) 

P*— /?c5 P’yc 

h$rn Qt 1 J 
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METHOD OF SOLUTION 


We have used the following unknowns for the expression of the flow field 


steady component 


unsteady component 


rz r 2 

f fts> 

r<«> , r/cs> 


v t?c<n , Tjfcn 

(fcrt , £sit > 


4r* , An i 

Aro », Arm , Ann, Arm 


These unknowns are determined by the following relations: 


B.C. on the volute: 
Continuity: 
Vorticity : 


Steady Component 
Eq, ( 24) 

Eqs. (32) ,(33) 
Eqs. (38) , (39) 


Unsteady Component 
Eqs. (25), (26) 

Eqs. (34). (35), (36), (37) 
Eqs. (41), (42), (43), (44) 


These equations include integrals related to the vortex distribution on the volute 
surface, which should be evaluated by some appropriate method. Equations (24) -(26) 
are integral equations for the vortex distributions on the volute surface and could 
be reduced to simultaneous linear equations by a singularity method. In the solu- 
tion of the vortex distributions the "Kutta condition" at the trailing edge should 
be applied. Strictly speaking the circulation around the volute fluctuates and a 
free vortex sheet is shed from the trailing edge of the volute. Since we are mainly 
interested in the forces on the impeller, we will neglect the effect of the free 
vortex sheet but apply the following conditions at the trailing edge. 


Steady part: 


II (so = 0 

645-) 

Unsteady parts: 


= r/(Sz) w (SO 

(46) 

~ to. S** -TdCSOh/CSt) 

can ) 


Now we can express all the relations as a set of simultaneous linear equations which 
can be solved numerically. The steady component may be solved independently of 
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unsteady component, and the result used in the analysis of the unsteady components 


UNSTEADY FORCES ON THE IMPELLER 


By considering the balance of the momentum of the fluid in the impeller, we can 
express the forces on the impeller as follows; 

Steady component 

x-zr ■ -i f - §c pi "'3 

+ ^ CuZ'- CU'-ZV') 1 ***' J C4?) 

and unsteady component 

(Xc~ Z&)CoSO>t t ( *1' 

-fc, $<*£'] + *)€’*”* 

^r[^CU'-^)CU^) ag' - f Ci ( U'~a tF' )CU'-jL9 r OdB'} 

■t i £ [§c z Cu'-itj') COS CB-Cot) Kale (U'-iU')coS(Q -O jCJ> r,ol&] 

-factf VrtK, e-fcrO r^ (r) e-^-f+P 3 r*rd& 

The total pressure is given by Equation (3) and the integrals can be evaluated 
analytically by using the expressions (28-31). 

CONCLUDING REMARKS 


The unsteady forces can eventually be expressed in the form of stiffness 
matrix, 

fX\ (MZ , Xs/z \fx*£e»sc*t 

M' \%/z , %lz )[ 

The time average of the force component in the direction of whirling motion is given 
by % (Y c - 5T s ) and the sign of this quantity determines whether or not^the fluid 
forces have destabilizing effects on the whirling motion. The sum Yt (x c +Y s ) gives 
the time average of the force component in the radial direction and thus the 
hydrodynamic stiffness. The ultimate goal of the present study is to examine these 
factors for realistic impeller-volute combinations. 
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